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Abstract 



(N 
O 

Cf) • A variational analysis of the pure SU(N) gauge theory in 3+1 dimensions at finite tem- 

perature is performed, extending the work of Kogan, Kovner and Milhano 0. A de-confining 
*"pH , phase transition is found at a temperature of 470 MeV, somewhat higher than lattice esti- 

mates |5J. This value is however rather sensitive, for reasons which are discussed. A more 
robust quantity is the ratio of the transition temperature to the lightest glueball mass in 
the model. This is 0.18, in agreement with the lattice estimate for SU(3) to two significant 
figures. Ways of further improving the calculation are discussed. 
PACS Numbers: 12.38.Aw, 12.38.Lg 
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1 Introduction 



In a recent paper P^, a variational method is used to study the deconfinement transition in the 
pure SU (N) gauge theory at finite temperature. The method mimics the Rayleigh-Ritz variational 
method in the Schrodinger formulation of quantum mechanics. There, the standard procedure is 
to take a physically motivated ansatz for the ground state wavefunction, parameterized by some 
free parameters, and to minimise the expectation value of the Hamiltonian with respect to those 
parameters. This provides an upper bound for the ground state (vacuum) energy. The method 
at finite temperature is analogous: in the canonical ensemble formulation of quantum statistical 
mechanics, one forms an ansatz for the density matrix, with free parameters, and minimises the 
expectation value of the Helmholtz free energy. This provides an upper bound for the free energy 
at a given temperature. 

In section|5J we begin with a precis of the approach followed in |J. The calculation generalises 
the variational analysis at zero temperature performed in [3]; an additional kernel H in the ansatz 
corresponds to the effect of thermal disorder in the system. This kernel is taken to be small, 
and only the leading order correction to the entropy of order H log H is considered. In this 
approximation a deconfining phase transition is found to occur at a temperature of 450 MeV. 

In section [31 we consider higher order corrections in H to the entropy. It is shown that 
(within approximations already present in the Kogan-Kovner model at zero temperature) one 
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can calculate the entropy to all orders in H in the high temperature phase. In this extended 
analysis, the transition temperature is shifted to 470 MeV, which is high compared to lattice 
estimates However, this value depends on a mean-field estimate of the critical coupling in 
a sigma model which arises in the analysis, and is therefore only approximate. The ratio of the 
transition temperature to the lightest glueball mass in the model, which is independent of this 
mean field estimate, is 0.18. This is in agreement with the lattice estimate for SU(3) to two 
significant figures. 

We conclude in section 01 by discussing our results, and suggesting further improvements. 



2 The order HlogH analysis 

The ansatz is constructed by considering density matrices which in the field basis have gaussian 
matrix elements, and where gauge invariance is explicitly imposed by projection onto the gauge- 
invariant sector of the Hilbert space. It reads 



p[A, A'} = J DU exp | - i AG' 1 A + K u G~ x X u - 2AHA' U X 



(1) 



where, under an SU(N) gauge transformation U, A — > A u and DU is the SU(N) group-invariant 
measure. In the above we employ a matrix notation, with e.g. 

AGHA = J dxdydz A1{x)G<${x - y)H%(y - z)A%(z) . (2) 

Here, indices i, j,k, . . . e {1, 2, 3} and a,b,c, . . . e {1, 2, . . . , iV 2 — 1} denote the spatial Lorentz 
components and colour components of the gauge field respectively. Explicitly, the gauge transfor- 
mations are 

AUx) - Af a (x) = S ab (x)A h i (x) + X}{x) , (3) 

with S ab = ±tr(T a U f T b U), A? = Hi-(T a WdiU), and ^ form an N x N Hermitian representation 

of SU(N): [^, ^] = if abc ^- with normalisation tr(r a r h ) = 2S ab . 

The kernels G _1 and H are arbitrary variational functions. To facilitate the calculation, they 
are restricted to be isotropic in colour and space indices. Furthermore, one splits the momenta 
into high and low modes with k ^ M and restricts the kernels to the one parameter momentum 
space forms 

The form for G _1 is motivated by the propagator for a massive scalar field, viz. (k 2 + M 2 ) 1 / 2 ; 
the form for H assumes that only the low modes are thermally excited at the temperatures of 
interest. 1 With the above restrictions on the kernels, only two variational parameters, M and H, 
remain. 

Before discussing the variational analysis at finite temperature, let us recall the analysis at zero 
temperature. The former will turn out to be a straightforward generalisation of the latter. At 
T = 0, H = and the analysis reduces to the minimisation of the energy, that is of the expectation 
value of the Hamiltonian, U = trTip/trp, where 

W=l[£ 2 + 5 2 ], (5) 

with Ef = 5/5A^ and Bf = e ljk (d 1 Al+gf abc A b J A c k /2). This is equivalent to the analysis originally 
performed by Kogan and Kovner in [3]. Firstly, one performs the Gaussian integrals over the gauge 

1 Non-zero H in Q corresponds to thermal disordering, since H = corresponds to a pure state. 
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fields A. This leaves integrals over the gauge transformations U, evaluated with respect to a sigma 
model 'action' which is both non-local and non-polynomial in U. To simplify the action, the 
gauge transformations U are split into parts dependent on high and low momentum modes, with 
k ^ M as above. The effect of integrating out the high modes is to effect a renormalisation group 
transformation: the coupling g 2 of the low mode sigma model is replaced by the renormalised 
coupling g 2 (M). So M acts as a UV cut-off for the low mode theory. Furthermore, the theory is 
asymptotically free @JE]. Provided M is sufficiently large (and g 2 (M) sufficiently small) one can 
then consider the low mode theory to leading order in g 2 (M) . The relevant Euclidean actions are 

S[A] = (a + £\ G- 1 U+£\ + jAG f_1 A , (6) 



r i U ] = TTTTT^^sutN) I d 3 x d i U\x)d i U[x) , (7) 



for the A fields and 

WW) xrsu{N) 

for the low mode U fields, where the trace is performed over SU(N) matrices U . The high modes 
do contribute to the energy at zero temperature. However, they do not yield any additional 
contribution at finite temperature, since H, which parameterizes the thermal disorder in the 
theory, is zero for k > M. 

Next consider this low mode sigma model as a statistical mechanical system at 'temperature' 
g (M). The system undergoes a phase transition with spontaneous symmetry breaking from a 
disordered state at small M (large g 2 (M)) to an ordered state at large M. Calculations in the 
disordered phase are performed in the mean field approximation: the U are treated as N 2 free 
fields obeying the unitarity constraint U* U = 1. In the ordered phase the sigma model is treated 
in leading order perturbation theory, writing U = e igv T I 2 and expanding the exponential. In 
the disordered phase, the energy 2 is minimised close to the phase transition with M ~ M c , 

N 2 M 4 

and g 2 (M c ) — tt 2 /N. In the ordered phase, one obtains 

rr N 2 M A 

so that the energy is indeed minimised at M ~ M c , on the disordered side of the sigma model 
phase transition. 

The extension to finite temperature was discussed in 1 . At finite temperatures, the energy 
minimisation argument is modified: one must consider the balance between energy U and entropy 
S, minimising the free energy F = U — TS. Since the parameter H corresponds to thermal 
disordering, one expects generically that S will vanish for vanishing H. In the SU (N) theory at 
moderate temperatures, the degrees of freedom correspond to glueballs. Since these are heavy, 
the excitations (disordering), and consequently the entropy, will be small. One can thus attempt 
to calculate the entropy as some expansion in the small parameter H. 

The leading order contribution in H is [J a term of the form H log H, multiplied by a coefficient 
which is an SU(N)l ® SU(N)r symmetric correlator of U fields. In the disordered (symmetric) 
phase of the sigma model, this expectation value vanishes. Furthermore, since the leading order 
contribution in H to the energy is positive definite, one finds that the free energy is minimised 
with H = at the minimum of the energy. 

Thus, in the disordered phase of the sigma model, the minimum of the free energy is at M ~ M c 
with, from ©, 

N 2 M? , . 



J Here and throughout, extensive quantities are written per unit volume. 
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In the ordered phase of the sigma model, the leading contribution to the entropy at small H is 



Then, from ©, 



N 2 M 3 

S =—^ Hl °SH- (11) 



N 2 M 4 N 2 M 3 , , . 



Minimising with respect to H and M, one finds that F is minimised in the sigma model disordered 
phase (with (U) = 0) from T = up to T = T c ~ 0.33M C , beyond which F is minimised with 
M in the ordered phase of the sigma model (with (U) ^ 0). Since U plays the same role as the 
Polyakov loop variable at finite temperature, this corresponds to a deconfinement phase transition 
in the pure SU (N) gauge theory. 

As a result of the minimisation procedure one finds that the dimensionless quantity H/M 
is equal to e — 1 . This raises the question of whether neglected terms of O(H), which have the 
same magnitude as the retained terms of 0(H log H), could considerably affect the calculation. It 
would, therefore, be desirable to extend the calculation to higher order in H. This we do in the 
remainder of this paper. 



3 Extended analysis 

We wish to extend the previous calculation to include terms beyond the leading order in the kernel 
corresponding to thermal disorder H. Although it is not at all clear a priori how one might do this 
in general, we shall, by providing an alternative procedure to that underlying the results of 013] 
outlined above, show that an extended analysis is indeed possible within the context and inherent 
approximations of the Kogan-Kovner model. It should be recalled that the stumbling block to 
any improvement is the calculation of the entropy, S — —trp log p. 

Let us try the following gambit. Instead of restricting ab initio the density matrix to the form 
(|T|). imagine that we take some arbitrary gauge-invariant density matrix ansatz depending on the 
A fields and integrated over the U fields. We allow this new ansatz (and whatever kernels it may 
contain) to remain arbitrary until we have no choice but to restrict it. Now we integrate out the 
A fields to obtain a partition function of U fields with respect to some action. 

Next we introduce a separation of momenta into high and low modes with k ;s M and integrate 
out the high mode U fields as before. This effects a renormalisation group transformation on the 
low modes, replacing the bare coupling g 2 — which is not arbitrary, since it is defined by the gauge 
transformations (JSJ — by the running coupling g (M). Now provided our ansatz is sufficiently 
close to the correct density matrix for SU(N), the theory will be asymptotically free. We are 
thus left with an action for the low modes which is again some complicated sigma model, with a 
renormalised coupling g 2 (M) which we expect to be small provided M is large and vice versa. 

Now consider this model as a statistical mechanical model at 'temperature' g 2 (M). We make 
the plausible assumption that this sigma model will, as M is varied, undergo a symmetry-breaking 
transition at 'temperature' g 2 (M c ) from a 'thermally disordered' (symmetric) phase at large 
g 2 (M c ) to an ordered phase at small g 2 (M c ). Further, it is clear — since the Polyakov loop (U) 
is zero in the former phase and non-zero in the latter — that this sigma model phase transition 
corresponds directly to the deconfinement transition in the SU(N) theory. 

This argument is quite general; on review, it is clear that our only assumptions are that the 
ansatz is sufficiently close to SU (N) and that the low mode sigma model undergoes a symmetry- 
breaking phase transition. In particular, let the ansatz, which is arbitrary and need not be 
Gaussian, be the correct density matrix for SU(N). The first assumption is certainly true. If the 
second assumption is also true, then we have constructed an exact argument that the deconfine- 
ment transition in SU(N) corresponds to the phase transition in the low mode sigma model. 

Thus, in order to study deconfinement in SU(N), our aim should be to model the physics of 
each sigma model phase as accurately as possible and calculate the transition scale M c . We then 
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calculate the free energy of SU(N) in each phase, including any possible contribution from the 
high modes, at temperature T and extract the minimal free energy. The deconfinement transition 
occurs at the temperature for which the free energies calculated in the ordered and disordered 
phases of the low mode sigma model coincide. 

Although we will take (JTJ as the ansatz for the density matrix, we shall keep the kernels G -1 
and H arbitrary until we have no choice but to restrict them. 

In the disordered phase no progress seems possible without restricting the arbitrary kernels. 
Following [J], we adopt the forms J2J as before and the analysis is identical. The Boltzmann factor 
is e~ Mfl / T in this case where M g is the lightest glueball mass, so we expand the small entropy to 
leading order and get zero as before. The resulting minimal free energy is thus independent of the 
temperature and we find 



F= — 



N 2 M* 
30tt 2 



(13) 



where M c ~ 1.33GeV is the sigma model transition scale predicted by the mean field calculation 
<>f EH- 

In the leading order perturbation theory approximation to the ordered phase of the sigma 
model, however, minimisation with respect to arbitrary kernels G~ 1 and H for both high and low 
modes is possible. Further, the analysis can, as desired, be carried out to all orders in the thermal 
disorder kernel H. 

In this approximation, the U matrices can be parameterised in the standard exponential form 
and expanded in the coupling g 



U = cxp<{ lj/y 



= 1 



W y + ■ 



(14) 



Hence at leading order one can take 



U~l, 

diU ~ igdup —■ 

Thus, the gauge transformations © reduce to 

At - At - d^ a 

and the Hamiltonian (JSJ reduces to 

n = \ [Ef + [e^Alf] 



(15) 
(16) 
(17) 



But these last two equations describe the theory U(1) N 1 : in the leading order of sigma model 
perturbation theory, the SU(N) Yang-Mills theory reduces to the U(1) N _1 free theory. Moreover, 
the density matrix JJJ becomes Gaussian again, because the gauge transformations are linear. One 
has 



[A, A']= J Dip exp | - i [AG' 1 A + (A' - d^G^iA' - dtp) - 2AH(A' - dip)] J 



(18) 



Now the theory of TV 2 — 1 U(l) free fields in 3 + 1 dimensions is completely tractable; the 
variational analysis for the U(l) theory (with Gaussian ansatz (|18fl ) was discussed in The free 
energy in momentum space in terms of the arbitrary kernels G" 1 and H is 



F 



N 2 



d 3 p 
(2tt)3 



G _1 (l + GH) +p 2 G(l - GH)- 1 



AT log 



GH 



(I - (GH) 2 ) 1 / 2 ~ (I - GH) 



- log 



1-(1-(GH Y) 
GH 



2)1/2 



1 - (1 - (GH) 2 ) 



2)1/2 



(l-(GH) 2 ) 1 / 2 -(l~GH), 



(19) 
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The kernels which minimise the free energy are 




1 — e t 

and the minimal value of the free energy at temperature T is 

„ N 2 -l 



(20) 



J™P 2 d P [| + ri g(i- e ^/ T ) 



(N 2 - 1)T 4 f°° , x 3 
1 ax 



3tt 2 

7r 2 (A^ 2 - i)r 4 

45 







(21) 



where the zero-point term has been discarded. All of this is of course consistent with the standard 
analysis of photon gases in statistical mechanics. 

The minimal free energy of SU (N) in the ordered phase of the sigma model at temperature T 
is, from l|21|) and dropping sub-leading contributions of 0(N°), 

n 2 N 2 T 4 , . 

F= . (22) 

45 v ' 

So we see that the free energy of SU (N) is minimised with M = M c in the disordered phase of 
the sigma model for temperatures from zero up to a temperature T c where 



which in turn implies 



N 2 M 4 ir 2 N 2 T 4 

30tt 2 45 ' V ' 



3\ 1/4 M r 



T c = [-) — - ~ 470MeV . (24) 



We note that the transition temperature is shifted by only a very small amount compared to the 
result T c ~ 450MeV obtained in pp. The calculation is improved in the sense that, in the high 
temperature phase of SU(N), which corresponds to the ordered phase of the sigma model, we 
have been able to extend the original analysis to include all orders of the thermal disorder kernel. 
This is desirable because at high T this kernel, which corresponds to the Boltzmann factor, is of 
order unity. The calculation is also improved in that the minimal kernels in the high T phase, 
approximated as free gluons, are the exact ones. If we had performed the calculation with the 
kernels (@J, we would not have been able to obtain the true minimum of the free energy. As 
in pQ, we find that the deconfinement phase transition is strongly first order with latent heat 
AE = ^iT c 4 . 

Finally, it will also be of interest to calculate the ratio of the transition temperature to the 
lightest glueball mass in the model, which is 2M C [7J- One obtains 

T i.&i/^O.W. (25) 



2M C 2tt k 2 

4 Discussion 

In this extended variational analysis, we have identified a phase transition (within the approx- 
imations made) at 470 MeV. This seems rather high in comparison with numerical simulations 
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performed on the lattice, which give around 280 MeV for 5/7(3) However, the estimate ob- 
tained for the transition temperature is only expected to be approximate since it is sensitive to 
the value of the critical scale M c , which is calculated in the mean field approximation of the sigma 
model. This sensitivity can be removed by computing the ratio of the transition temperature to 
the lightest glueball mass in the model (|25|l . One then obtains a value of 0.18, which agrees with 
the lattice result for SU (3) to two significant figures. 

There are other reasons why we expect the calculation to be only approximate. The most 
important point to be aware of is that in the original zero temperature analysis, the SU(N) gauge 
theory was hived into two parts (the high and low modes) for the purpose of tractability. The 
former corresponds approximately to the perturbative gauge theory, which is well understood (and 
which we have treated in the zeroth order) and the latter to the low energy sector, which is less 
well understood and is treated in the mean field approximation. In considering the theory at finite 
temperature, the phase transition corresponds to a jump between the two sectors. So in doing the 
analysis, we are really asking the question: at what temperature does the free U(1) N _1 gauge 
theory become thermodynamically more favourable than the low energy theory calculated in the 
mean field approximation? 

Whilst this is a perfectly sensible question, to which we have obtained a sensible answer, one 
must ask whether this means anything for the full SU (N) gauge theory. One is interpolating 
between a low energy theory, which is already only approximate, and a high energy theory which 
is only correct in the ultra-violet limit. This is, to say the least, rather crude. However, we believe 
that the principle of the method is rather powerful, in that there is scope to improve upon the 
calculation. The simplest way in which this can be done is to include perturbative corrections to 
the free energy coming from the high modes. The finite-temperature corrections should be added 
to HJ . In contrast, only the zero-temperature corrections should be added to (|13|l . since there 
are no thermal contributions in this sector by assumption. 

The second way in which improvements can be made is to improve corrections to the low 
mode sector. We believe that the crudest approximation here is in taking the leading order of 
perturbation theory in the ordered phase of the sigma model. Clearly this is not appropriate close 
to the sigma model phase transition, which corresponds also to the SU(N) phase transition in 
this model. So if there are large shifts in the sigma model behaviour near the transition, we would 
expect the transition temperature to be significantly shifted. An improved treatment of the sigma 
model near the phase transition necessarily calls for a higher order or non-perturbative calculation 
to be performed. But then one has to calculate the entropy for a non-free theory. Such a task is 
beyond our present calculational abilities. However, we are encouraged by the fact that the all 
orders in H result for the free energy in the ordered phase obtained in this paper differs only very 
slightly near the phase transition from the one obtained in [I] where only the leading HlogH 
contribution was taken into account. This is, of course, why the transition temperature is not 
significantly shifted by the improved analysis. Now a non-perturbative calculation of the entropy 
to order HlogH does seem to be possible, and is currently under way [Hj- 
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